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SCALING  CALCULATIONS  FOR  A  RELATIVISTIC  GYROTRON 


I.  Introduction: 

The  gyrotron  oscillator  is  under  development  as  a  source  of  high  power 
millimeter-wave  radiation.  Research  in  the  U.  S.  to  date  has  emphasized  high 

average  power  devices  based  on  thermionic  cathode  technology.  Output  powers  up  to  a 
megawatt  have  been  achieved  by  such  devices  which  are  based  on  high  quality  moderate 
voltage  (<  100  keV)  electron  beams  with  currents  of  <  50  Amps.  High  efficiencies 
~  50%  have  beeen  demonstrated.  Future  directed  energy  applications  of  high  power 
millimeter-wave  radiation  may  require  sources  with  peak  power  in  the  gigawatt 
range.  It  is  therefore  of  interest  to  investigate  the  peak  power  potential  of  a 
relativistic  gyrotron  with  a  mega-el ectronvolt,  mi Iti-ki lomampere  electron  beam. 
Such  an  electron  beam  can  be  generated  for  short  pulse  lengths  (<  100  nsec)  using  a 
pulseline  accelerator  with  a  field  emission  (cold)  cathode.  An  important  objective 
of  this  work  is  to  obtain  design  parameters  for  a  cold  cathode  gyrotron  experiment 
based  on  a  600  kV,  6kA,  55  nsec  Febetron  pulser.  A  second  objective  is  to 
investigate  scaling  to  higher  voltage  to  increase  output  power. 

Section  2  of  this  report  describes  scaling  parameters  and  design  constraints 
for  a  relativistic  gyrotron.  The  results  of  calculations  are  given  in  Section  3. 
The  non-linear,  si ow-time-scal e  equations  of  motion  used  for  these  calculations  are 

derived  in  the  Appendix. 
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II.  Gyrotron  Scaling  Theory  and  Design  Constraints 

This  section  describes  scaling  parameters  and  constraints  for  a  relativistic 
gyrotron  with  a  megavolt,  multi-ki loampere  electron  beam.  The  configuration 
analyzed  corresponds  to  a  cylindrical  resonator  and  a  thin  annular  electron  beam 
with  the  beam  radius  chosen  to  coincide  with  a  maximum  of  the  resonator  electric 
field. 

The  gyrotron  is  a  particular  case  of  the  Cyclotron  Resonance  Maser  (CRM).  The 
resonance  condition  is 

w-k  v  -sn/y-0  (1) 

n  11 

wnere  w  is  the  rf  frequency,  fi  is  the  non-relati vi stic  cyclotron  frequency,  y  is  the 
relativistic  mass  ratio,  v„  is  the  axial  electron  velocity,  k  is  the  axial 
propagation  number,  and  s  is  the  harmonic  number.  In  the  gyrotron  the  electron  beam 
interacts  most  effectively  with  a  TE  type  mode  which  is  close  to  cut-off  in  the 

resonator.  Thus,  the  gyrotron  corresponds  to  the  case  k  v  <<  to  -  so/y  .  The 

present  calculations  consider  only  the  fundamental  harmonic  interaction  (s=l)  as 
this  yields  the  highest  power  and  is  least  subject  to  mode  competition. 

The  basic  gyrotron  interaction  is  characterized  by  five  parameters:  the  e-beam 
relativistic  factor  y  ,  the  e-beam  pitch  ratio  a  =  v  /v  ;  the  number  of  cyclotron 
orbits  an  electron  makes  traversing  the  resonator,  Nc;  the  applied  axial  magnetic 
field,  B0;  and  tne  amplitude  of  the  rf  electric  field  at  the  e-beam  position,  Eb. 
Gyrotron  designs  for  high  power  are  subject  to  constraints  on  the  minimum  resonator 
0,  the  maximum  electric  fields  at  the  cavity  wall,  the  space-charge  limit  on  beam 
current  and  requirements  for  stability  of  single  mode  operation.  To  optimize  the 
gyrotron  design  for  a  given  e-beam  energy,  characteri zed  by  y ,  and  pitch  ratio  a  , 
it  is  convenient  to  use  the  number  of  cyclotron  orbits  Nc  as  the  scaling 

parameter.  The  applied  magnetic  field,  B0,  and  the  electric  field  amplitude  at  the 

e-beam  are  then  varied  to  optimize  the  efficiency. 


In  terms  of  Nc,  the  interaction  length  is  given  by 


L  =  —  N  =  (N  +  1)  X  8,rt 
(I»  JO  C  C  #0 


(2) 


where  Eq.  (5)  below  has  been  used  and 
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so 


2 

a 


(3) 


is  the  normalized  axial  velocity,  u>c  is  the  relativistic  cyclotron  frequency 
(u>c=  and  x  =  — is  the  free-space  wavelength.  The  rf  fields  in  a  cylindrical 

cavity  can  be  expressed  in  terms  of  the  axial  profile  function  h(z).  The  profile 
function  depends  on  several  factors  including  the  interaction  itself;  however,  a 
reasonable  approximation  for  scaling  calculations  is  to  assume  a  sinusoidal  profile: 


h(z)  =sink||z 


(4) 


where  kB  =  tt/L.  The  interaction  efficiency  is  optimized  when  the  cyclotron 
frequency  is  detuned  from  the  resonance  frequency  by  approximately  the  interaction 
bandwidth,  i.e.. 


«c 

“  '  uc  *  W  (5) 
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Thus  a  good  estimate  for  the  optimum  magnetic  field  is  given  by 


B  = 
o 


fef  1  +  1/N 


(6) 


which  was  used  in  these  calculations.  This  leaves  the  electric  field  amplitude  Eb 
as  the  only  optimization  parameter. 


Estimates  for  the  electronic  efficiency  were  obtained  by  integrating  a  set  of 
equations  of  motion  for  an  electron  in  the  cavity  and  averaging  over  the  initial 
phase  of  the  electron  entering  the  cavity.  A  tractable  set  of  non-linear  equations 
based  on  a  slow-time-scale  formulation  developed  previously  was  used.  For  this 
application,  the  equations  derived  in  reference  1  for  a  TE  circular  waveguide  mode 
have  been  generalized  to  include  the  interaction  with  the  rf  magnetic  f’eld 
component,  which  is  important  at  relativistic  e-beam  energies.  A  derivation  of  the 
equations  of  motion  is  given  in  the  Appendix.  Electron  beam  self-field  effects  and 
thermal  effects  are  neglected. 

As  shown  in  the  Appendix,  the  slow-time-scale  equations  of  motion  for  a 
gyrotron  interacting  at  the  sth  harmonic  with  a  TE  circular  waveguide  mode  which  is 
close  to  cutoff  can  be  written  in  the  form: 
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In  Equations  (7)  u^  and  a  are  slow-time-scale  variables  for  the  magnitude  and  phase 

of  the  transverse  momentum,  and  uz  is  the  axial  momentum.  c  is  the  axial 

2  2  2 

coordinate  and  y  is  the  relativistic  factor  (y  =  1  +  u^  +  uz).  E$  is  the 
effective  electric  field  at  the  beam  and  h  is  the  normalized  axial  profile  function 
of  the  electromagnetic  mode  ff=  |f  |  h(z)l.  J$_^  is  a  regular  Bessel  function. 


-  A  A*  A* 
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a)  is  the  normalized  wave  frequency  and  nzQ  is  the  nonrelati vistic  cyclotron 
frequency  due  to  the  applied  axial  magnetic  field.  Equations  (7)  use  the  following 
normalization  scheme: 


(i)  the  radial  coordinate  is  normalized  to  r  (r  =  r/r  ),  and  the  axial 

wo  wo  * 

coordinate  is  normalized  to  r*w0/xmn  (z  =  z  xmn/rwoh  where  x(I)n  is  the  nth  zero  of 
dJm/dx; 

Oir 

(ii)  frequencies  are  normalized  to  c/(xmnrw0)  («  =  _^2_) ; 

mn 


(iii)  velocities  are  normalized  to  c  and  momenta  are  normalized  to  mQc 
(  V  =  v/c,  u  =  p/moc); 

2 

me  x 

(iv)  electric  and  magnetic  fields  are  normalized  to  -i — t—  and 

'  '  wo 

m  c  x  iir  |  i  r 

Tel  "  r^-’  respectively  (E  =  -*--U  £,  8  Unnormalized 

i  I  -wo  mn  oc  mn 


quantities  are  expressed  in  mks  units  unless  otherwise  stated. 

In  obtaining  Equations  (7),  the  radial  scale  length  rWQ  is  assumed  to  be  equal 
to  the  waveguide  radius  and  the  axial  dependence  of  the  waveguide  radius  has  been 
neglected  except  in  calculating  the  axial  profile  function  f. 

Equations  (7)  are  independent  of  the  wave  frequency.  Moreover,  geometric 
factors  associated  with  a  particular  waveguide  mode  and  the  position  of  the  electron 
beam  are  contained  in  the  factor  E$. 

“  i  me 

For  a  circularly  polarized  mode  with  azimuthal  dependence  e 
(and  time  dependence  e1<ot  )  E$  is  given  by 


.(c)  Jm-s 


(kmn  SJ 
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and  by 
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w 

for  a  mode  with  dependence  e+im0.  in  Equations  (8)  and  (9)  r^  is  the  normalized 
waveguide  radius;  Rq  ,  =Q  are  the  polar  coordinates  of  the  electron  orbit  guiding 
center  orbit  and 


r(c)  = 
mn 


2  2  1/2 
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The  profile  function  f  has  dimensions  of  voltage  and  is  normalized  according  to 

f  -  Mr  The  transverse  wave  number  k  =  x  /r 
2  mn  mn  w 

me 
o 

For  a  linearly  polarized  mode  with  azimuthal  dependence  cos  m  o 
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where 
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The  output  power  is  related  to  the  cavity  fields  according  to 


P  = 
(0 


u)W 
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where  the  stored  electromagnetic  energy 
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(13) 


(14) 
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and  0  is  the  diffraction  Q  of  the  cavity.  Substituting  Eq.  (14)  into  Eq.  (13)  and 
transforming  to  normalized  parameters  leads  to 


.  3  m  ‘  - 

P  =  —  —  (—  1  —IT  I 
to  2  to  e  0  1  max' 


2  2 
fdc  j  h ( z )  | 


where  c  is  the  speed  of  light  and  uQ=  4ttx10" 


Eq.  (15)  shows  that  the  output 


power  is  independent  of  frequency.  The  electronic  efficiency  is  calculated  by 
integrating  Eqs.  (7)  through  the  interaction  region  to  obtain  the  change  in  electron 
energy  and  averaging  with  respect  to  the  initial  momentum  phase: 


n  =  r0dA0[Y  (cF.  A01  -  T0]  /  (Y  -1) 


For  a  linearly  polarized  mode  the  efficiency  must  also  be  averaged  with  respect  to 
the  guiding  center  angle  =Q.  The  electron  beam  power  or  current  are  found  from 
the  power  balance  equation: 


Ph=  IV-  - 
b  o  o  n 


The  threshold  current  for  oscillation  is  given  by  the  small  signal  limit 

P  (E  ) 

[th=  lim  TTT T)  (18) 

th  E$+  0  V'V 

Calculations  of  gyrotron  efficiency  generally  assume  a  circularly  polarized  rf 
mode  because  in  this  case  the  interaction  efficiency  is  the  same  for  all  beam 
electrons.  When  the  rf  mode  is  linearly  polarized  the  rf  fields  have  an  azimuthal 
dependence  given  by  sinme  or  cosme  and  consequently  the  electronic  efficiency  varies 


with  azimuthal  angle.  As  shown  in  Section  3,  the  result  is  a  reduction  in  the 
optimum  efficiency  averaged  over  the  beam  by  a  factor  of  approximately  0.7  compared 
to  a  circularly  polarized  node.  Linear  mode  polarization  is  nevertheless  of 
particular  interest  for  relativistic  gyrotrons  as  this  polarization  occurs  when 
axial  slots  are  used  to  control  mode  competition. 

In  these  calculations,  a  value  of  twice  the  diffraction  limited  0  was 
assumed.  The  diffraction  limited  0  is  given  by 


%  -  4’<V  1)2 


In  gyrotrons  the  space-charge  limited  electron  beam  current  places  a 
significant  restriction  on  the  achievable  output  power.  An  estimate  for  a  thin 
annular  gyro-beam  has  been  derived  by  Ganguly  and  Chu^.  Their  result  for  the  space- 
charge  limited  current  is 


iscl^<a)  = 


rfir? ' {1  "*<> )1/3]' 

in  (  w/  G) 


where  rw  is  the  cavity  wall  radius  and  r^.  is  the  e-beam  radius. 


III.  Results  of  Calculations 

This  section  presents  scaling  calculations  for  the  relativistic  gyrotron.  To 
show  the  efficiency  potential  of  the  relativistic  regime  and  to  calibrate  the 
present  calculations.  Figure  1  compares  efficiencies  calculated  for  70  keV  and  1  MeV 
electron  beams  with  a  =  v  /v^  =  1.5  interacting  with  a  circularly  polarized  TE 
mode.  One  notes  that  the  peak  efficiency  of  30%  at  1  MeV  is  only  slightly  less  than 
the  peak  efficiency  of  34%  at  70  keV.  However,  there  is  an  order  of  magnitude 
difference  in  the  optimum  number  of  cyclotron  orbits:  peak  efficiency  at  1  MeV 
occurs  for  Nc  =  3-5  whereas  peak  efficiency  at  70  keV  occurs  for  Nc  =  15-30.  This 
result  is  consistent  with  the  gyrotron  efficiency  scaling  relation 


PvD 

c  o 


derived  by  Bratman  et  al  .  The  efficiency  results  for  70  keV  are  in  good  agreement 
with  the  results  of  Chu  et  al4  The  latter  are  also  shown  in  Figure  1. 

Although  the  intrinsic  efficiency  has  only  a  gradual  decrease  with  increase  in 

beam  energy  other  factors  limit  efficiency  at  relativistic  energies:  The  efficiency 

is  a  sensitive  function  of  the  ratio  v^v^  since  as  shown  by  Bratman  et  al3 
2 

n  a  8ln  Due  to  lack  of  experience  in  the  design  of  high  voltage  beam  formation 
systems  -  particularly  cold  cathode  systems  -  lower  values  of  v./v  can  be 
expected,  at  least  initially.  As  discussed  above,  the  use  of  a  linearly  polarized 
rf  mode  also  reduces  efficiency.  Figure  2  shows  efficiency  as  a  function  of  Nc 
for  a  600  keV  e-beam  for  a  =  1.0  and  1.5  and  for  linear  and  circular 

polarization.  For  a  linearly  polarized  mode  and  a  51  1.0  -1.5  the  peak  efficiency 
is  15%  -  22%  compared  to  21%  -  31%  for  a  circularly  polarized  mode.  Similar  results 
for  a  1  MeV  beam  are  shown  in  Figure  3.  Effects  due  to  electron  beam  thermal  and 
self-field  effects  have  not  been  included.  These  effects  may  be  significant  in  a 
high  power,  relativistic  gyrotron  and  are  under  investigation. 


Figure  4  shows  the  optimum  electric  field  vs.  Nc  both  at  the  electron  beam 
position  and  at  the  cavity  wall  for  a  linearly  polarized  mode  and  a  600  keV  beam 
with  a  =  1.5.  The  electric  field  corresponds  to  the  case  of  a  TE14  mode  with  the 
electron  beam  placed  at  the  fourth  E-field  maximum  (rG/rw  =  0.86).  The  curves  shown 
in  Figure  4  indicate  that  the  peak  cavity  field  at  the  wall  is  reduced  an  order  of 
magnitude  from  the  field  at  the  beam.  Thus  peak  fields  at  the  beam  of 
>  1  MV/cm,  which  are  needed  for  high  power  and  peak  efficiency,  can  be  achieved  with 
low  probability  of  field  induced  emission  from  the  cavity  wall. 

Figure  5  shows  the  electron  heam  current  corresponding  to  optimum  efficiency  as 
a  function  of  Nc  for  a  600  keV  beam.  The  cavity  0  is  assumed  to  be  twice  the 
diffraction  limited  0  and  the  beam  is  placed  on  the  fourth  E-field  maximum.  Results 
are  shown  for  the  TE14  mode  with  a  =  1.0  and  1.5,  and  for  the  TE15  mode  with  a  = 
1.0.  Linear  polarization  is  assumed  in  all  cases.  The  estimated  space-charge 
limiting  currents  for  these  configurations  are:  2.9  kA  for  the  TE15  mode  with  a  = 
1.0,  3.5  kA  for  the  TE14  mode  with  a  =  1.5  and  7.4  kA  for  the  TE}4  mode  with  a  = 
1.0. 

Output  power  scaling  at  optimum  efficiency  for  a  600  keV  beam  is  shown  in 
Figure  6.  The  choices  for  0  and  beam  position  are  the  same  as  in  Figure  5.  The 
maximum  rf  power  is  limited  by  the  dc  space  charge  effect  for  the  TE14  mode,  a  = 
1.5  and  TE^  mode,  a  *  1.0  cases.  The  maximum  power  for  the  TEj4,  a  3  1.0  case  is 

limited  by  the  available  e-beam  current  of  6  kA  corresponding  to  a  Febetron 

pulser.  Figure  7  compares  output  power  scaling  for  1  MeV  and  0.6  MeV  e-beams. 

The  curves  correspond  to  a  TE^4  linearly  polarized  mode  with  a  *  1.5,  the  beam  on 

the  4th  E-field  peak,  and  0  =  20n.  The  1  MeV  curve  reaches  a  peak  power  of  about  1 

GW.  The  dashed  curve  indicates  that  higher  powers  are  possible  if  partial  space- 
charge  neutralization  proves  feasible. 


The  present  analysis  has  emphasized  the  high  power  regime.  The  interaction 
region  is  very  short  involving  only  a  few  cyclotron  orbits,  and  non-linear  effects 
are  strong.  In  this  regime  the  output  power  is  ultimately  limited  by  the  condition 
that  the  operating  e-beam  current  he  greater  than  the  oscillation  threshold 
current.  Figure  8  compares  the  threshold  and  optimum  efficiency  currents  for  a  600 
keV  beam  with  a  =  1.5  placed  on  the  fourth  peak  of  a  TE^  mode.  The  optimum 
efficiency  currents  are  shown  for  both  linear  and  circular  polarization;  the  former 
being  10%  -  15%  larger  than  the  latter.  Roth  polarizations  have  the  same  threshold 
current.  Except  at  short  interaction  lengths  of  less  than  3  cyclotron  orbits,  the 
optimum  efficiency  current  is  approximately  twice  the  threshold  current.  The  same 
behavior  is  found  for  a  1  MeV  beam. 

In  the  case  of  short  pulse  operation  it  is  necessary  to  consider  the  output 
power  risetime.  The  10%  -  90%  power  risetime  is  approximately 

r  -  —  [  \ - 1  ]  1  (22) 

“  thr 

where  I  is  the  beam  current  and  Ithr  is  the  threshold  current. 

This  expression  does  not  include  the  effect  of  finite  e-beam  risetime. 

Figure  9  shows  power  risetime  scaling  vs.  Nc  for  0.6  MeV  and  1  MeV  beams  with  ct  = 
1.5,  0  =  20^.  The  minimum  risetime  for  a  35  GHz  oscillator  is  shown  to  be  about  one 


nanosecond. 


IV.  Conclusions 

The  results  of  this  study  indicate  that  the  relativistic  gyrotron  has  potential 
for  achieving  high  efficiency  (15-30%)  using  relativistic  (y~2-3)  electron  beams. 
Optimum  efficiency  occurs  for  short  interaction  lengths  characterized  by  4-8 
cyclotron  periods.  Such  short  interaction  lengths  lead  to  the  possibility  of  very 
high  peak  power  generation  using  mult i-ki loampere  beams:  100-300  MW  for  a  600  keV 
beam  and  ~  1  GW  with  a  1  MeV  beam.  These  power  levels  apply  to  the  case  of  single 
mode  operation  in  a  linearly  polarized  TE14  mode. 

The  optimum  interaction  efficiency  occurs  for  electric  fields  ~  1  MV/cm  at  the 
electron  beam.  Since  the  gyrotron  is  based  on  a  fast-wave  interaction,  such  fields 
at  the  beam  can  be  achieved  without  excessive  electric  fields  at  the  wall.  For 
example,  for  a  600  keV  beam  and  a  TE^4  mode  the  peak  field  at  the  wall  is  ~  100 
kV/cm. 

The  output  power  rise  time  has  been  estimated  to  be  a  few  nanoseconds  for  a  low 
Q  oscillator  at  35  GHz.  Thus  the  gyrotron  interaction  can  be  readily  investigated 
using  a  20-100  nsec  pulseline  accelerator. 

The  electron  beam  dc  space-charge  effect  appears  to  be  a  significant  factor 
limiting  gyrotron  power.  The  possibility  of  space-charge  neutralization  via  a 
background  plasma  should  be  investigated. 
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Append! x 


Sinql e-Particle  Slow-Time-scale  Equations  of  notion 


The  starting  point  is  the  Lorentz  force  equation 


Ht  (yy)  +  ^  (y  x  K 


n  (2  +  v  x  §)  =  a 


where  §0  is  an  applied  static  magnetic  field  which  is  predominantly  axial 


§  =  z  +  3  r  ,  8  >>  B  , 

o  oz  or  oz  or 


=  M 


is  the  charge  to  rest  mass  ratio  for  an  electron,  and  t,  §  are  the  electromagnetic 
fields.  It  is  convenient  to  introduce  the  normalized  momentum  variable 


■*  «► 
U  =  yV 


in  terms  of  which  the  relativistic  factor  is  given  by 


Y  -  [1  +  U2/C2]l/2 


To  obtain  a  si ow-time-scal e  formulation  the  transverse  momentum  is  expressed  in  the 


u  +  i  u  ■  i  u  exp  [i  (  n  t  +  <fr)  ] 
x  y  c 


(6) 


where 


n 


pB  (z*o) 
oz 


(7) 


and  r  *  t  -  t  where  t  is  the  time  the  electron  enters  the  interaction  region. 

o  o 

Using  the  transformation  (6)  in  Eq .  (1)  leads  to  the  following  equations: 


if 


(a  -  ia  ( 

x  y 


a  +  i  a  )e  T  J 

x  y 


u  u 
z  z 


JLLB 

B  3z  (8, 


1) 


f  (a  -ia 
x  y 


)  i(nx+4>) 


+(a  +  ia 
x 


(8.2) 


2 

11  3B  Ut 

11  a  A  —  —  —  -  -  -  ■ 

z  z  2  B  3z  y 


(8.3) 


Equations  8. 1-8. 3  account  for  electron  interactions  with  all  rf-field  components  and 
allow  for  a  tapered  applied  magnetic  field. 


2 .  Transverse  Electromagnetic  modes  in  a  waveguide 

In  the  present  theory,  the  electron  beam  is  assumed  to  interact  with  a  single 
TE  or  TM  vacuum  waveguide  mode  with  time  dependence  elwt.  The  transverse  field 
components  are  given  by: 
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?t=  Re  (f(z)  e  (x,y , z)  eiwt} 


(9.1) 


lt  =  Re  {g(z)  &  (x,v,z)  e^*1} 


(9.2) 


where  e  and  £  are  transverse  vector  mode  functions .  These  vector  functions  are 
related  to  scalar  functions  according  to 


e  *  z  x  V 


h  -  -  Vc  i, 


for  TE  modes,  and 


+ 

e 


-  Vt  * 


+ 

z  x  V  t|» 


for  TM  modes,  The  scalar  functions  satisfy  the  equation 


(7*  +  kL2)  «  -  0 


and  the  boundary  condition 


=  o 
3n 


on  the  waveguide  boundary  for  TE  modes,  or 


(10.1) 


(10.2) 


(11.1) 


(11.2) 


(12) 


(13.1) 
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at  the  waveguide  boundary  for  TM  modes.  In  Eq .  (13.1)  3/3n  indicates  the  normal 
derivitive.  When  integrated  over  the  waveguide  cross  section  the  vector  functions 
satisfy  the  normalization  and  orthogonality  conditions 


/S(z)dxdy  Vj  *  6i.1  (14> 

-► 

(5..  ■  0  if  i*j  ,  5  =  1  if  i  *  j)  where  v  denotes  either  e  or  ti  and  S(z) 

denotes  the  z-dependent  waveguide  cross  section.  The  axial  field  components  are 
given  by 


B  -  Re  {-i  —  f(z)  (x,y)  eiljC} 
z  1  w 


(15) 


for  a  TE  mode,  and  by 


E  = 
z 


2  2 


Re  {-i  -  k  g(z)  if)  (x,y)elutj 
1  u)  l 


(16) 


for  a  TM  mode.  The  axial  profile  functions  f(z)  and  g(z)  are  related  according  to 


g(z) 


i  djf(z) 
to  dz 


(17) 


for  a  TE  mode,  and 


for  a  TM  mode.  Equations  (15-18)  are  obtained  by  substitution  into  the  Maxwell  curl 


equation 


..  f! 


3.  Equations  of  motion  for  an  electron  interacting  with  a  waveguide  mode 


Substituting  Equation  (9)  -  (16)  into  (8)  leads  to 


ut=  Re  { -n  (f-vzg)  (L^e-1  ^  +  L  j>  e1 


(&T+$  )  . 


“zut  l_^z 
~TT  Bz  3z 


l!  (■^■r(f.,zg)  L_*e1,aT**’) 


.2  V^_f* ,  e  i.t,  .  a  +  h. 

ui  1  y 


;  21*  ,  at*  e-1'0'**’  ♦  L_*  e*  <»**♦> 


^ _ 1_  i^z. 

2y  Bz  3Z 


for  TE  modes,  and 


at»  Re(-1n(f-v2g)  tnT4*'-i_,e1 (n,4» 'le1”'! 


Vt  1  sRz 


Re  {  -  |  (f-vzg)  (L+*e 


■i  (Or+$)  +  L 


ei(nT^))(?iut^ 


li7 

+—  (21.1) 

r 


uz  =  Re  {-  -in.  g  rvt(L+te-i(nT+,>)-  Lj*1  (flT+* }) 

2  2  2 

c  S  it  i  3R7  ut 

*  ^  -r  •  i  -  *57  sr  f-  <*'•*> 

for  TM  modes,  where  L+  and  L  _  are  defined  by: 

M  =  ei9  (Jr  +7fe  <**•" 

L  _tj»  =  e“i0  (  ~  "r  |©  )  ♦* 

3.1  Circular  waveguide  mode. 

For  a  circular  waveguide,  the  projection  of  an  electron  orbit  on  the  cross 
sectional  plane  of  the  resonator  forms  an  annulus  with  average  radius  equal  to  the 
electron  orbit  guiding  center  radius  RQ,  and  with  thickness  equal  to  twice  the 
Larmor  radius 


rL  5  \/az 


(23) 


The  case  of  an  axis  encircling  beam  (RQ=0)  is  a  special  case  of  the  analysis. 

The  scalar  function  for  a  circular  waveguide  mode  with  circular  polarization  is 


=  C  J  (k  r)e 
mn  mn  mn 


time 


(24) 


where  for  a  TE  mode: 


(26) 


C  =  [[it  (x 
mn  mn 


-  m‘)] 


i/c 


Jn,(xmn>  1 


-l 


xmn  is  the  ntn  zero  of  Jm  (x)  (prime  denotes  differentiation),  Jm  is  a  Bessel 
function,  and  rw  is  the  waveguide  wall  radius.  For  a  TM  mode: 


mn 


mn 


Cmn  =  [  /ir  zmn  Jm  (znJ  ] 


(27) 

(28) 


and  zm  is  the  nth  zero  of  Jm(z).  for  definitness  consider  the  mode  with 

polarization  e^ut  ”  IT,e^.  Substitution  of  Fq.  (?4)  into  (22.1)  and  (22.2)  gives 


L  (J)  3  Cm„  kmn  l(k™r)e 

+  mn  mn  m-1  mn 


-i (m-l)0 


L  =  C  k  J  .  (k  r)  e 
-  mn  mn  m+1  mn 


-i (m+l)0 


(29.1) 

(29.2) 


For  the  case  of  non-zero  guiding  center  radius,  the  addition  theorem5  is  used  to 
express  Bessel  functions  of  the  radial  coordinate  in  terms  of  Bessel  functions  with 
guiding  center  and  Larmor  radius  arguments;  i.e.. 


J„(kr) 


+ime  v 
e  =  > 


wv 


Jk(*rL)et,kU-(!lTt‘>i,kw',i  1  (30) 


where  =Q  is  guiding  center  position  angle. 

Substitution  of  Equations  (24),  (29),  and  (30)  into  (20)  or  (.21),  and  considering 
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Che  interaction  with  the  sth  harmonic  of  the  applied  magnetic  field  leads  to  the 


following  slow  time-scale  equations  of  motion:  For  a  TE  mode: 


3 J  (k  r  ) 

u  -  Re  {-n  (f-v  g)  C  J  (k  R)  —2.  mn 
t  1  z*  mn  m-s  mn  6  3  rT 

L 


i  [  (ui-sft  )t  +  ait  -  s$  -  (m-s)H  ]  i  ,  UzUt  1  8Bz 
0  0  }  + ~TT  r 


2y  B  3z 
'  z 


(31.1) 
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(31  .2) 


u  .  3J  (kmnrL) 
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z  1  y  mn  m-s  3r_ 


1 1  n  ^ 

i[(ai  -  SOT  +  ait  -  sA  -  (m-s)  =0]  ,  1  t  1  _ z  . 

e  0  '  2  y  Bz  3z 


(31.3) 


and  for  a  TM  mode: 
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In  Eqs.  (7)  in  Section  2,  above  the  phase  angle  $  is  replaced  by 


A=(u)-sfi)T+wT  -  s$-(m-s)=  . 

o  o 


(33) 
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